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ABSTRACT 

Snap roll is a phenomenon which has plagued the sub- 
marine designer since the advent of the high speed submarine. 
A description of snap roll and its possible causes is pre- 
sented in this thesis. A computer program is developed 
which simulates a submarine in surge, sway, yaw and roll. 

It is verified that snap roll is very sensitive to the 
metacentric height. Additionally, it is seen that reducing 
the sail size does not have as large an effect on reducing 
the roll angle as anticipated. Rudder sequencing and speed 
reduction are two other methods examined. The results of 
the four investigations are presented in graphic form. The 
conclusions include a justification for the use of the 
simulation computer model as a design tool. 
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NOMENCLATURE 



S YMBOL 



CB 

CG 

D 




J 

L 



m 



N 



N 



i 



P 



DEFINITION 
Center of buoyancy. 

Center of gravity. 

Drag force. 

Moment of inertia about the i axis. 

Product of inertia about jk axis. 

Hydrodynamic moment about x axis 
(rolling moment). 

Non-dimensional coefficient used in 
representing K as a function of i . 

Overall ship length. 

Lift force. 

Ship’s mass. 

Hydrodynamic moment component about 
z axis (yawing moment). 

Non-dimensional coefficient used in 
representing N as a function of i . 

Angular velocity about the x axis. 
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Angular acceleration about the x axis. 

Angular velocity about the y axis. 

Angular acceleration about the y axis. 

Angular velocity about the z axis. 

Angular acceleration about the z axis. 

Velocity of the origin of the body axes 
relative to the fluid. 

U velocity component in the x direction. 

Acceleration component in the x 
direction. 

Initial steady value of U in the x 
direction. 

U velocity component in the y direction. 

Acceleration component in the y 
direction. 

U velocity component in the z direction. 

Acceleration component in the z 
direction. 

The x coordinate of the center of buoyancy. 
The x coordinate of the center of gravity. 
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The x coordinate of the center of 
pressure of the sail. 



Hydrodynamic force along x axis 
(axial force). 



Non-dimensional coefficient used in 
representing X as a function of i. 



The y coordinate of the center of 
buoyancy. 



The y coordinate of the center of 
gravity. 



The y coordinate of the center of 
pressure of the sail. 



Hydrodynamic force along y axis 
(lateral force). 



Non-dimensional coefficient used in 
representing Y as a function of i . 



The z coordinate of the center of 
buoyancy . 



The z coordinate of the center of 
gravity. 



The z coordinate of the center of 
pressure of the sail. 



Angle of drift. 



Bowplane or sailplane deflection. 



Rudder deflection. 



Sternplane deflection. 

Ratio u /u. 
o 

Angle of pitch. 

Angle of yaw. 

Angle of roll. 

Constants used to represent the 
propeller thrust in the axial 
equation. 



CHAPTER I - INTRODUCTION 



1.1 BACKGROUND 

Since the advent of the high speed submarine in the 
late fifties, the submarine designer has been confronted 
with new design problems in the area of ship control. One 
of these problems, which is still with us two decades later, 
is snap roll . Part of the reason for the perpetuation of 
this problem has been a general lack of adequate data and 
information concerning this phenomenon. Consequently, 
appropriate design criteria could not be developed and the 
problem of snap roll has continued. Naval architects and 
ship designers need adequate alternatives and trade-offs to 
provide good designs. This is the motivation for this 
thesis . 

A computer model will be developed in Chapter II which 
simulates a submarine in roll, yaw, surge and sway. This 
model can be used by naval architects as a design tool in 
performing trade-off studies or as a method of establishing 
design criteria for a particular design. The alternatives 
for reducing snap roll, which are investigated in Chapter III 
using this simulation model, were based on the immediate 
needs of the design community and the time and monetary 
constraints normally placed on this type of effort. It was 
felt that these alternatives should not increase the com- 
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plexity of present designs or alter present design 
practices, in order that they may be put to immediate use. 
Chapter IV presents the conclusions and recommendations of 
this thesis. In the appendix will be found a listing of 
the computer model, a flow chart and a sample output. 



1.2 WHAT IS SNAP ROLL ? 

An excellent description of snap roll was presented by 
Griffin, et al ( 1 ) ^ : 



"At 20 sec after full rudder, the initial 
transients have died out, and the lateral 
states begin to follow the trim values 
which correspond to the decreasing forward 
speed. Eventually, the forward velocity 
reaches its trim value, and r, q, v, w 
and are fixed at their respective trim 
values. For r, q, v, and w, the difference 
between the plateau value and the later trim 
value is not large; however, for roll, it 
can differ by more than 13 deg. This is 
snap roll . " 



Simply stated, snap roll is the maximum roll amplitude 
that results from a turning maneuver. The large roll angle, 

^ Numbers in parentheses indicate references listed 
at the end of this thesis 
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which develops early in the turn, constitutes the "snap." 

The large peak value of the roll angle is the submarine 
tending towards its steady state turning value for the 
initially high forward speed. This usually occurs about 
20 seconds into the turn. As the submarine enters the 
turn and during the initial phase of the turn, there is a 
large side slip velocity. This side slip velocity probably 
accounts for a major portion of the maximum roll angle. 

The subsequent reduction of the large roll angle is pri- 
marily due to speed loss. The high drag force developed 
during the turn maneuver dramatically reduces the forward 
velocity component. The steady state forward speed can be 
expected to be 1/4 to 1/3 of the initial speed. 

There is still conjecture amongst hydrodynamicists as 
to the theory and physics involved in explaining this 
phenomenon. It is generally agreed that snap roll pri- 
marily depends on the turning rate, the amount of rudder 
deflection, the fairwater (sail) size and location, and the 
ship speed. A large turning rate and/or rudder deflection 
produces a large lateral velocity. Consequently, this large 
v acting on the sail produces a large roll moment. Addi- 
tionally, it has been proposed that the starting lift 
exceeds the steady lift due to the dominating effect of the 
trailing vortex sheet downwash in the steady case (2). 

This may be due to the starting vortex formed and shed by 
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the sail. However, there is no substantiating evidence for 
this . 
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CHAPTER II - THE SIMULATION MODEL AND 
COMPUTER PROGRAM 

II. 1 THE FOUR DEGREE OF FREEDOM MODEL 

There exists at least three computer programs which 
simulate a submarine in a maneuver. However, none of these 
programs are available for general use, primarily due to 
proprietary considerations. Therefore, in order to inves- 
tigate snap roll in submarines, it was necessary to develop 
another simulation model. Unlike the other models, however, 
this model involves the use of four degrees of freedom vice 
six degrees of freedom equations. The four coupled equations 
considered here are the X, Y, N and K equations. The justi- 
fication for the use of only four equations of motion is 
that in a turn, depth excursions are small and the Z and M 
equations may be decoupled. The depth change which occurs 
during a turn, called "squat”, normally can be easily 
remedied with a small deflection of the diving planes and/or 
sternplanes. These small excursions have a negligble effect 
on snap roll. This is evidenced by the fact that a four 
degree of freedom simulation program provides valid results 
for a turning maneuver. This will be seen from the results 
presented later in this chapter and in Chapter III. The 
advantages afforded by a four degree of freedom program, 
in addition to providing valid results, are: 
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• Reduced computation time. 

• The estimation or measurement of fewer 
hydrodynamic coefficients. 

Of course, such a model, necessarily restricts one's con- 
sideration to only maneuvers involving surge, sway, roll 
and yaw. 

II. 2 THE COORDINATE SYSTEM AND EQUATIONS 

The coordinate system chosen for use in this model was 
based on reference ( 3 )* This reference gives the standard 
equations of motion for submarines used by the U.S. Navy. 

The coordinate system is illustrated in figure 2-1. It 
should be noted that the origin is centered at the center of 
gravity. A better choice may have been to use the center of 
buoyancy; however, to avoid a rewriting of already valid 
equations, the reference ( 3 ) system was used. The center of 
buoyancy of the basic hull for a submarine design is 
necessarily the geometric center of the hull. Therefore, a 
transverse velocity (or acceleration) would not produce a 
roll moment, caused by the basic hull form and this would 
reduce the number of necessary calculations. 

The equations used in the development of this simu- 
lation model were taken directly from reference (3) and are 
listed in Appendix B. There are no second order or higher 
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order acceleration terms involved in these equations. This 
is due to the lack of any significant interaction of the 
viscous and inertia forces. This allows the use of potential 
theory which gives adequate values for the hydrodynamic 
forces (4). Acceleration and velocity interaction terms 
are also assumed negligibly small for the same reason stated 
above. It should be noted that the equations remain non- 
linear even without these second order acceleration terms, 
due to the higher order velocity terms. 

The equations are applied to a submarine in a turning 
maneuver under the following assumptions and restrictions: 

1. The submarine is neutrally buoyant, i e. 

7/ - B = 0. 

2. The submarine is at least three hull 
diameters below the free surface, i.e. 
fully submerged . 

3- The z component of the velocity is zero. 

4. There is no initial list, i.e. ^ = 0. 

5. There is zero trim, i.e. 9 = q = q = 0. 

6. The origin is at the center of gravity, 
i-e- X G Z y G = Z G = °- 

7- Only the moments of inertia along the 
principal axes are non-zero, i.e. 

I xy = I y7i ~ I xz = 0. This is partly a 
consequence of #6. 
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8. 



The coefficients and equations used are 
in non-dimensional form. 

9. The control surface deflections are 
initially zero, with only the rudder 
angle being non-zero for t>0. 

II. 3 THE SOLUTION OF THE EQUATIONS OF MOTION 

The simplest method of solving for the roll angle 
induced by a turning maneuver is to solve for u(t) from 
the non-linear X equation; then, to use this u in solving 
the linear Y, N and K equations. This will provide com- 
paratively accurate values for the roll angle for relatively 
small deflections and for small forward speed losses. How- 
ever, we are interested in high speed and tight turns, 
where snap roll is most prominent. The vehicle motions in 
this type of turning maneuver are outside the valid linear 
range. Therefore, the non-linear equations of motion are 
used. 

The general method used for solving the four non- 
linear equations is an iterative technique. The specific 
technique employed in this thesis was adapted from the 
method used in reference (5)- The four equations are each 
written in terms of the accelerations u, v, r and p. They 
take the form*. 
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X: 



(m‘- X^) u - y G r = f x (u, v, r, p, £ r ) 

Y: (m’- Y.p v - (m’z G + Y^ 1) p + 

(m’x G - Y^ 1) = f 2 (u, v, r, p, S r ) 

N: (i; - H! ) r - (i; x + N- ) p + 

(m’x G /l 2 - N^/l ) v - (m’y G A 2 ) u = 

(u, v, r, p, S r ) 

K: (i; - K! ) p - (i; z + Kp r - 

(m’z G /l 2 + K^/l) v = f 4 (u, v, r, p, S r ). 

The prime indicates a non-dimensional value. 

Expressions for the accelerations u, v, r and 
p are derived by employing matrix methods. Based on the 
initial conditions u Q , v , etc. , starting values for u, v 
r and p are obtained for t =0. Values for u, v, r and p 
are calculated from these starting values employing the 
Taylor expansions, 

u(t + At) - u(t) + At*u(t) 
v(t + At) - v(t) + &t-v(t) 
r(t + At) - r(t) + At-r(t) 

^(t +At) - ^(t) + &t*p(t) + ( (At) 2 )- p( t) . 



These new velocity values are then used to update the pre 
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viously calculated accelerations. This indicates the 
iterative nature of this integration technique. The value 
of S r and t are updated after each iteration. The equations 
are integrated over the time span specified. The size of At 
is arbitrary; however, the smaller the At, the greater the 
accuracy obtained from this technique. The accuracy desired 
should be tempered by the fact that small At's require 
rather large amounts of computer time. It was found that 
At’s on the order of one half to one second provided very 
satisfactory results, while not using excessive amounts of 
computer time. 

II. 4 COMPUTER SIMULATION MODEL VALIDITY 

A large number of computer runs were executed for 
various submarines, under various conditions. It was in- 
tended to prove the validity of the developed computerized 
simulation model using hydrodynamic coefficients taken from 
towing tank model tests. The results from one of these 
computer runs is presented in figure 2-2. As can be seen, 
in general, the model accurately predicts the submarine’s 
response in roll; however, the predicted snap roll angle is 
too large. This can be attributed to a number of causes. 

One is the possibility of scale effects, since the co- 
efficients used in the program are not the full scale 
hydrodynamic coefficients. Secondly, the metacentric height 
used was that specified in publications. The actual GM 
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Figure 



during the full scale trials could vary slightly from this 
value and, as will be seen in Chapter III, the maximum roll 
angle is very sensitive to small changes in GM. However, for 
other submarines and sequences, the predicted maximum roll 
angles were less than the full scale values and therefore, 
this somewhat indicates that the difficulty in predicting 
the snap roll angle does not lie entirely in the model or 
technique used. It will also be noted that the model 
results increasingly deviate from the full-scale trials as 
t increases. This can be attributed to two causes. First, 
the full-scale trials are performed with a spiraling 
maneuver. In order to save computer time and since we were 
only interested in the maximum roll angle, the computer 
simulation runs were only performed for 180 0 turns. There- 
fore, the rudder angle was taken off earlier for the model 
and the roll angle decreased more rapidly for large t. 
Secondly, at slower speeds (large t), the cross coupling 
effect of pitch and heave on roll is more pronounced than at 
higher speeds. Since we are using only four degree of 
freedom equations, the effect of pitch and heave at low 
speeds is lost. 
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CHAPTER III - METHODS FOR REDUCING SNAPROLL 



111.1 GENERAL 

There are a number of techniques and methods which have 
been proposed to reduce snap roll in submarines (1). Four of 
these alternatives will be investigated here. As was in- 
dicated in the introduction, the choice of these alter- 
natives was based on the needs of the design community and 
time constraints. The four methods chosen can be categorized 
under one of two headings. The first concerns changes in the 
naval architectural characteristics of the design. Increasing 
GM and reducing sail size are classed under this category. 

The second category involves those alternatives which make 
use of an automatic ship control system. Rudder sequencing 
and speed reduction fall under this classification. Each 
method of snap roll reduction will be investigated indi- 
vidually in the following sections. 

111. 2 METACENTRIC HEIGHT 

The metacentric height is directly related to the 
righting moment, K^, by the relation, 

K<p - -W • GM 

where W is the submerged displacement. This term is some- 
what disguised in the governing equations in Appendix A. 

In the K-equation, 
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K<p — Zg • B 

where B is the buoyancy, and equal to the submerged dis- 
placement at neutral buoyancy, and Zg is the metacentric 
height. It can be shown that, as a first approximation, the 
maximum roll angle is inversely proportional to the meta- 
centric height (1). Therefore, it would be advantageous to 
increase the metacentric height in order to reduce the 
effect of snap roll. However, this could prove to be a 
difficult design task, since it would involve adding weight 
or rearranging weights in an already weight-limited design. 

The metacentric height for several submarines was 
increased by 0.25> 0.50, 0*75> 1.0 and 1-5 feet. A non- 
dimensional plot of the results are shown in figure 3-1. 

It should be noted that the amount of rudder used and u 

o 

were different for each vehicle. This indicates that the 
percent decrease in the maximum roll angle is essentially 
independent of u o and the maximum ordered rudder angle for 
a given GM change. The curve is asymptotic for greater than 
100$ changes in metacentric height. Therefore, increases in 
the metacentric height of greater than about 100$ will yield 
no greater than approximately an 82$ decrease in the maximum 
roll angle. 

As an example, a 30$ increase in GM yields a 33% de- 
crease in . For a typical submarine, this translates 

into a 10° decrease in the maximum roll angle for an approx- 
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imately 4 inch increase in the metacentric height. It can be 
seen that large decreases in the snap roll angle can be 
realized for relatively small changes in GM. This is 
evident in figure 3-2. This is a plot for a typical sub- 
marine which shows the decrease in the maximum roll angle 
for increasing GM. An average 5* 75 ° decrease in is 

IT13.X 

gained for each 3 inch increase in GM. A point of minor 
interest shown in this plot is that the maximum roll angle 
occurs earlier in the turn as GM increases. 

III. 3 RUDDER SEQUENCING 

Rudder sequencing means to defer the full rudder angle 
ordered, until some predetermined speed loss has been reached. 
It is expected that by limiting the initial rudder angle, 
the maximum roll angle will be reduced. This is based on 
the fact that the roll angle is directly proportional to the 
rudder deflection. Rudder sequencing is envisioned as part 
of an automatic ship control system; however, the method of 
rudder sequencing is simple enough to be employed as a 
manual process. One disadvantage of this method of snap roll 
reduction is that we might expect the transfer to increase. 
However, such an increase is normally no more than one half 
to one ship's length. The procedure employed in this thesis 
to investigate this roll reduction method involves reducing 
the maximum ordered rudder angle by 5» 10 , 15 and 20 degrees. 
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Effect of Reducing GM on a Typical Submarine 
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Figure 
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The predetermined speed loss criteria is to maintain the 
reduced rudder deflection until the speed in the turn has 
been reduced to one third that of the initial speed into the 
turn. This is depicted in figure 3-3- The selection of 
this speed criteria was based on the fact that the steady 
speed in a fully developed turn is approximately one third 
of the initial speed. The rudder deflection rate was held 
constant for all submarines and all sequences at 4 deg/sec. 

The results are depicted in figure 1-^. As can be seen 
from the plot, this method of snap roll reduction is effec- 
tive. However, relatively large rudder angle reductions are 
needed to produce small changes in . For example, a 10% 

reduction in the maximum roll angle requires a 23 f reduction 
in the initial rudder angle. The use of rudder sequencing 
for a typical submarine is shown in figure '3-5. It can be 

seen that a greater decrease in *P per degree of rudder is 
o max J a 

gained for the 20° rudder sequence than for the 30° rudder 
sequence. It should be noted that some of the figure 3-U 
data points were based on the same submarine at different 
initial speeds, i.e. 24 and 30 knots. Therefore, rudder 
sequencing seems to be independent of speed, in the sense 
that the same benefit is received f or a particular rudder 
sequence despite the initial speed. For example, if a 20 
degree rudder sequence reduces snap roll by 15 % for u q = 20 kts, 

then the same sequence will also reduce snap roll by 15% for 

a u = 30 kts. 
o J 
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Rudder Sequence for a Typical Submarine 
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Figure 3-5 



III. 4 SAIL SIZE 



A submarine sail produces a large roll moment which is 
caused by the sideslip velocity. The magnitude of the roll 
moment varies widely from submarine to submarine, primarily 
due to the extensive variation in sail location and size. 

The sail dimensions are primarily derived from considerations 
for the quantity of equipment and the dimensions of the 
equipment to be placed in the sail. With the exception of 
the fairwater, which encloses the sail, the naval architect 
pays little attention to the hydrodynamic considerations of 
sail design. It was felt that if the sail size could be 
reduced, that a corresponding reduction in snap roll would 
be realized. 

An investigation was performed to determine if the 
chord or the span had the greater impact on snap roll. 

Figure 3-6 shows the results. Obviously, a change in the 
span will have a greater impact on roll, as evidenced by 
the larger slopes for the three coefficients considered. 

This result was anticipated since the moment arm associated 
with the spanwise location of the center of pressure will 
have a greater influence on roll than the moment arm 
associated with the chordwise location. Additionally, a 
reduction in span will produce a corresponding reduction in 
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the lift due to the lower aspect ratio. A decrease in the 
chord would obviously increase the aspect ratio. Conse- 
quently, the span was reduced by 1, 5* 10* 1 5 * 20 and 30 
percent for a number of submarines. The results from these 
computer runs will be presented; however, first the method 
used to predict the impact on the hydrodynamic coefficients 
of reducing the span will be presented. 

The approximate method for determining the impact of 
the sail on the hydrodynamic coefficients was partially 
derived from Abkowitz's presentation for a lifting foil in 
reference (4). This approximate method involves determining 
the Y force, N moment and K moment at the sail as a result of 
velocity disturbances v, r, and p. These coefficients are 
reduced correspondingly for a reduction in the sail size. 

The new sail coefficients are added to the original co- 
efficients, less the original sail contribution. Then the 
equations of motion are solved for T. Figure 3-7 defines 
the coordinate system, the variables involved in the 
derivation, and sign conventions. 

For most purposes, the sail x and y velocity components 
are the same as the submarines’s x and y velocity components, 
i . e . 



u - u 



s 



v 



= v . 

s 
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The subscript "s" denotes the local condition at the 
sail. The angle of attack produced at the sail by the 
transverse velocity is 

/3 = tan -1 v s = tan -1 - . 

Consequently, we can calculate the Y force and N moment 
produced at the sail as 

- (L cos + D si n/h) 




where L and D are the lift and drag forces on the sail, 
respectively. In terms of their respective non-dimensional 
coefficients, the lift and drag are 



L = C 



L ’ 



D => C 



D ' 



i A s u s 2 = C L ■ if A s (u s 2 + v s 2) 
i/ A s (u s 2 t v s 2 ) 



where A g is the projected sail area and j> is the fluid 
density. 

From airfoil theory (6), the slope of the lift co- 
efficient curve versus the angle of attack is shown to be 
approximated by 
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<) C L = 2 TT 



3/3 = 1 + 2/AR 



AR is the "reflected aspect ratio." This is equal to twice 
the mean span divided by the mean chord, which corrects for 
the proper pressure distribution on the sail. Using the 
above relations, it can be shown that the Y force on the 
sail due to a lateral velocity, v, is 



For the effect of an angular velocity, v, where 

v = x . r, we find 
s s 




and the N moment due to v is 




Similarly, for the K moment 




z 



s 



Y due to r 
s 



(Y ) .x . r 
v v' s s 



N s due to r * (Y y ) g • x g 2 • r 



and 



K s due to r = (Y y ) s . x s . z s . r 
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Taking the derivative of these coefficients with respect to 

r, 



(Y ) = (Y ) -x 

v r s vs s 

(N ) = (Y ) • x 2 

V r ' s ' v'S S 

(K ) = (Y ) • x • z 

x r s vs s s 

A similar process for the contribution of an angular velocity 

with v = -z • p, yields 
s s 

<Vs = “'Vs ■ z s 

<Vs 1 -<Vs ' x s ' 2 s 

<Vs = '<Vs • z s 2 • 



A similar derivation provides the hydrodynamic coefficients 
for the sail with respect to the accelerations v, r and p. 

The key coefficient is, of course, (Y^.) . Based on the added 
mass calculation for a flat plate of dimensions s/2 and c, 

(Y v ) s is 



< Y v>s 
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This form of (Y- ) is used with a correction factor de- 

V s 

veloped by Abkowitz (1977) and finally becomes 



(Y- ) 
vs 




0. 54 

(1 + s/c + c/s) 



where s is the mean "reflected span" and c is the mean 

chord. The mean "reflected span" is directly related to 

the "reflected aspect ratio", in that, it corrects for the 

proper pressure distribution on the sail. 

Since x and z are the distances from the center of 
s s 

pressure to the center of gravity along their respective 
axes, the following relations (7) proved beneficial, 



(CP)- = 0.25 * 3 
(cp) - = (4/3) 7T- s . 

(CP)- is the distance to the center of pressure chordwise 
location from the sail's leading edge. Similarly, (CP)^ 
is the distance from the root chord to the center of 
pressure . 

Employing the above described procedure and equations, 
new coefficients were developed for reduction in sail size 
for a number of submarines. The results presented in 
figure 3-8 are for two different submarines, with what are 
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Change in the Snap Roll Angle due to a Reduction in Span 
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Figure 3-3 



termed large and small spans. "Large" and "small" relate to 
the relative size of the sails when compared to similar 
submarine designs. It can be seen that a greater benefit 
is realized by the larger sail. The maximum reduction 
applied to the span was restricted to 30 %. Any larger 
reduction in span would not be feasible due to physical 
restrictions, such as minimum periscope height requirements. 

The decrease in the maximum roll angle through sail 
size reduction was not as large as was expected. It is 
believed that this may be explained by the fact that the 
sail contributes to the roll damping at the same time that 
it is acting as a roll producer. Therefore, by reducing 
the sail size, the roll damping is reduced, as well as the 
roll production. Fortunately, as indicated by figure 3-8, 
roll damping appears to decrease less rapidly than roll 
production for a given sail size reduction. 

III. 5 SPEED REDUCTION 

It has been well established that the forward velocity, 

u, has a strong influence on the roll angle when in a turn. 

It can be shown that the maximum roll angle is proportional 

to the square of the velocity. Therefore, it is expected 

that small speed reductions will give significant reductions 

in . This is an important point, since if a submarine 

max 
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is designed for a maximum speed in excess of 30 kts, then 
the operating community would like to be able to take full 
advantage of this speed in various maneuvers , as well as 
in a straight line motion. The greatest difficulty with 
snap roll occurs at a submarine’s top speed. It would be 
attractive, from an operational point of view, to be able 
to enter a turn at almost maximum speed. With this in 
mind, several submarines were simulated to enter turns at 
from one to four knots below their maximum speed. The 
results are presented in figure 3-9- As was expected, 
small changes in the velocity produce relatively large 
reductions in snap roll. For example, a 10%> reduction in 
u, decreases the maximum roll angle by almost 20%. This 
method of reducing the effect of snap roll is equally 
feasible for use in an automatic control system or as a 
manual method of operation. 
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CHAPTER IV - CONCLUSIONS AND RECOMMENDATIONS 



IV. 1 CONCLUS IONS AND RECOMME'NDATIONS 

From this study, it must he concluded that we can, with 
sufficient engineering accuracy, predict full-scale sub- 
marine dynamic performance, at least in roll, from model 
data. This is accomplished through the use of a com outer 
simulation model. This is of importance to the submarine 
designer for three reasons. First, this prediction tech- 
nique provides the ship designer with an invaluable design 
tool. The stability characteristics of a multi-million 
dollar submarine design can be accurately predicted prior to 
making a major commitment to the procurement effort. The 
ship designer now can proceed with confidence on a design, 
at least from a dynamic stability standpoint. Secondly, the 
computer model can be used by the submarine designer to 
perturbate and iterate on a particular baseline design. 
Consequently, such a technique could be employed to optimize 
a design with respect to its stability criteria. Finally, 
this method of predicting stability performance can be 
employed by the designer to establish design criteria. 
Without a foreknowledge of a submarine's performance in a 
particular maneuver, the ship designer can test past 
similar designs to establish some pattern of performance 
and thereby derive a stability performance criteria. 
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It is apparent from the results of Chapter III that a 
small additional amount of GM in any submarine design will 
significantly reduce snap roll. Bearing this in mind, it 
would seem appropriate to suggest that the U.S. Navy revise 
its present design criteria for GM. The present criteria 
requires a minimum metacentric height, which does not work 
well for all designs from a stability standpoint. It is 
proposed that the revised design criteria for GM be based on 
a maximum acceptable snap roll angle. This would be a more 
flexible design criteria which could be tailored to indi- 
vidual designs. More importantly, it would couple a per- 
formance measure to the naval architectural characteristics. 
It would also accomodate for the variations in size of the 
roll moment producers which are encountered from design to 
design, i.e. the larger the roll moment producers, the 
larger the required GM to minimize the effect of snap roll. 

The results of the other methods considered in Chapter 
III suggest obvious measures to be taken in an effort to 
reduce snap roll; these are to reduce speed, reduce initial 
rudder angle and reduce sail size. As was indicated in 
Section III. 4, the reduction of sail size was not as 
effective in reducing snap roll as was anticipated. The 
possible reason for this was also suggested in Section III. 4. 
However, a more thorough study in this area is certainly 
warranted. This is particularly true, in view of the fact 
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that an approximate method was used for predicting the 
appropriate hydrodynamic coefficients for the contribution 
of the sail to snap roll. A more precise method of pre- 
diction, coupled with a careful consideration of the total 
effect of the sail and sailplanes on roll stability, could 
lead to a better insight into the causes and minimization 
of snap roll. Other recommendations for future investi- 
gations are: 

• The use of differential sailplanes and/or 
sternplanes as roll control fins. 

• The use of unbalanced rudders. 

• The effect of increasing the automation of 
submarine control systems on stability 
characteristics . 

Throe areas of improvement for the computer simulation 
model are proposed. The first is the extension of the 
present model to use the six degrees of freedom equations of 
motion. This should not only improve the accuracy of the 
results obtained from the present four degree of freedom 
program, but it would also allow for the investigation of 
more complex maneuvers and dynamic reactions, e.g. squat 
with trim plus roll. Secondly, the model should be modified 
to permit for more ease in its use by submarine designers. 
The less complex the input requirements and the more lucid 
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the output, the wider use the computer simulation model will 
experience within the design community. Lastly, it is 
recommended that the present program and future versions be 
optimized with respect to computer time. With the in- 
creasing expense and the time-sharing requirements of com- 
puter facilities, wider use of this type of computer 
simulation model will be realized by an efficient use of 
available computer time. 
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A . 1 FLOW CHART 
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coefficients for 
differential 
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accelerations 
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Read I . C . ' s 



<! 




* heading 



Print 





51 





52 




Initialize 
and 

^iimensiona.liz) 
variables . 



Test 

or appropriate 
xthrus t/drj 
;oef£- 



Calculate 
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of differential 
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Return 
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A. 2 COMPUTER PROGRAM 
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A. 4 INPUT/OUTPUT VARIABLES 



VARIABLE 



AI, BI, Cl 



Al, A2 



B 



COURSE 



DELB 



DELGI' 



DSLRO 



^ ND = non 



MEANING 



FORMAT DIMENSIONS 



INPUT/ 

OUTPUT 



Set of con- 
stants repre- 
senting the 
propellar 
thrust in the 
X-equation 



Limits used 
for selecting 
the proper 
propeller 
thrust 



Ship's 

buoyancy 



Amount of 
ordered course 
change 

Initial sail- DIO. 4 deg I 

plane 

deflection 

Amount the DIO. 4 ft I 

original GM 
is to be 
changed 



Rudder de- D9-3 deg 

flection at 
time T 



dimensional 



DIO. 4 ND I 



DIO. 4 ND I 



DIO. 4 lbs I 

DIO. 4 deg I, 0 
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DELS 

DZLT 

ID 

INDEX 

IX, IY, IZ 
IXY , IXZ, IYZ 

KCOEFF 

L 

Ivl 

NCOEFF 



Initial DIO. 4 deg I 

sternplanc 

deflection 



Time increment DIO. 4 sec I, 0 

used in 

iteration 



40 character A1 ND I, 0 

alpha- 
numeric 
heading 



If >0, read 12 ND I 

new submarine 
coefficients . 

If — 0 , run 
same sub- 
marine for new 
initial 
conditions . 



Moments of DIO. 4 ND I 

inertia 



K-equation DIO 

coefficients 

Ship’s over- DIO 

all length 

Ship’s mass DIO 

N-equation DIO 

coefficients 



4 ND I 

4 ft I 

4 ND I 

4 ND I 
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P. Q» R 


Initial 
angular 
velocity 
about the x , 
y and z axes, 
respectively 


DIO. 4 


deg/sec 


I 


PDT , QDT, RDT 


Initial 
angular 
accelerations 
about the x, 
y and z axes, 
respectively 


D10.4 


deg/sec^ 


I 


PD TO , RDTO 


Angular 
accelerations 
about the x 
and z axes , 
respectively 


D9-3 


deg/sec^ 


0 


PHI 


Initial angle 
of roll 


DIO. 4 


deg 


I 


PHIO 


Roll angle at 
time T 


D9-3 


deg 


0 


PSI 


Initial angle 
of yaw 


DIO. 4 


deg 


I 


PS 10 


Yaw angle at 
time T 


D9-3 


deg 


0 


RHO 


Sea water 
density 


DIO. 4 


slugs/f 


I 


RO 


Angular 

velocity about 
the z axis at 
time T 


D9-3 


deg/sec 


0 
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RRATE 


Average 
rudder rate 


DIO. 4 


deg/sec 


I, o 


rudai:t 


Maximum 

ordered 

rudder 

deflection 


DIO. 4 


dog 


1,0 


RUD1 


Value of 
reduced 
rudder 
deflection 
associated 
with sub- 
routine ALTRUD 


DIO. 4 


deg 


I, o 


T 


Present time 


D9-3 


sec 


0 


THETA 


Initial angle 
of pitch 


DIO. 4 


deg 


I 


TLAG 


Time lag of 
of the control 
system 


DIO. 4 


sec 


I 


u, V, w 


Initial 
forward , 
lateral and 
vertical 
velocities 


DIO. 4 


kts 


I 


UDT , VDT , ?'DT 


Initial 
forward , 
lateral and 
vertical 
accelerations 


DIO. 4 


kts/hr 


I 
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UDTO, VDTO 


Forward and 
lateral 
acceleration , 
respectively 
at time T 


D93 


kts/hr 


0 


(JO 




Initial 

forward 

velocity 


D9-3 


kts 


0 


U00, VO 




Forward and 
lateral 
velocities , 
respectively 
at time T 


D9-3 


kts 


0 


WT 




Ship's weight 


DIO. 4 


lbs 


I 


X3, YB, 


ZB 


The x, y, z 
position of 
the center of 
buoyancy 


DIO. 4 


ft 


I 


XCOEFF 




X-equation 

coefficients 


DIO. 4 


ND 


I 


XG, YG, 


ZG 


The x, y, z 
position of 
the center of 
gravity 


DIO. 4 


ft 


I 


XX, YY 




The moving x 
and y coordi- 
nates of the 
submarine' s 
origin v/ith 
respect to a 
fixed coordi- 
nate system, 
whose origin 
is the point 


D9-3 


yds 


0 



where the turn 
begins 
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YCOEFF 



Y-equation DIO. 4 

coefficients 



ND I 
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APPENDIX B 



B.l AXIAL FORCE 



m 



[u - 



2 2 

vr + wq - x r (q + r ) t y P (pq - r) + z„(pr + q) 



|i 4 [fqq ' q2 * X rr ' r2 +X rp' r p] 



u 3 h'* 



* x vr' vr + X wq ' W P 



'] 




+ £ f 2 [x uu 'u 2 + X vv -v 2 + I w 'w 2 



£r u 2 [ X $rSv ' fir2 * X * s is SsZ - X Sb Sb ’ 5b2 



+ 1 t 
2 A 



a 



2 2l 

.u+b.uu +c.u 
1 l c l c J 



- (W - B) sin 0 






X vv 7 ’ v2 + x wv, 7 '" 2 * x SrSr 7 '^ V 



♦ (7- 1) 
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B.2 LATERAL FORCE 



m v - wp 




2 p 

+ ur - y G (r + p) + z Q (qr - p) + x Q (qp + r) 



Y-'r + Y p 'p t Y p|p jp|p| t Y pq 'pq ♦ Y qr 'qr 



Y- ’v + Y • vq + Y 'wp + Y 
v vq H wp 1 wr 



’ wrj 



V Ur + Y p 'up + Y |r( ^ r ’u|r|S r 



Y v|r |*xl (v2 - w2)2 



2 I I r ij 



Y*’u 2 + Y v 'uv t Y y j y | 'v \(v Z + w^) 2 | 



2 . ...2sh 



Y 'vw t Yx ' vfi&v 
vw Or 



B) cos 0 sin ¥ 



Y r? ’ur(7 - 1) 



Y V 7 ' UV + Y v|v| 7 ' v l (y2 * + Y Sr7'^r u2 J 



(7 - 1) 
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I 



roK> 



B. 3 ROLLING MOMENT 



V + - y qr - (r + pq^xz + (r2 - q2)I yz + (pr - q)I xy 



+ m 



y n (w - uq + vp) - z p ( v - wp + ur) 



'U 



K i'P + K r' r * K qr' qr + K »o' Pq * K dId|'- d Ip| 



pq ^ T "p|p| 



+ |jt 



4 



K ’ ud f K 'ur + K • ' v 
p * ' r v 



4 

n K ’ vq t K ’ wp + K ' wr 
V vq H wp ^ wr 



♦! i 



K*'u 2 K„' uv + K„,..rv|(v 2 + w 2 )^' 



v 



4 I 

V I V 



♦Si 3 



K ’ vw + Kr 'u^6r 
w i or 



] 



f (y„W - y B B) cos 9 cosf’- ( z n W - z R B) cos 9 sin f 



G B 



t || 3 k» 7 ’u 2 (7 - i) 
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B . 4 YAWING MOMENT 



+ (I y - I X )PQ - (4 + r P) I yz + (q 2 " P 2)I X y + (rq “ P^x 



+ rn 



x G (v - wp + ur) - y G (u - vr + wq)| = 



3= 



+ N p'P + Npq'PI + V' qr + N r|r|' r l r| 

lrf f v 



I /["P' 



v + N ' wr + N ' wp 4- N ' vq 
wr v/p - 1 - vq ^ 



up 4 N^ ' ur + 



N |r|<£r’ u * r ^ r + N | v |r’ l (y2 4 w 2 > T l^] 



+ | / 3 [ N ^’ u2 + N v’ UV + N v|v«' v |(y2 + w2) *'l 



+ Il 3 



N ' vw + Nr ' u^&r 
vw or 



+ ( Xq W - x B B) cos 9 sin<p + (y Q W - y fi B) sin 



9 



4 |i rf Ur ( ^" 1} 



+ fl 



’k,'Pv + N v|v |^v|(v 2 + w 2 P| + Nj^'^u 2 



(? - 1) 
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